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We present a mass map reconstructed from weak gravitational lensing shear measurements over 139 deg2 
from the Dark Energy Survey (DES) Science Verification data. The mass map prohes both luminous and dark 
matter, thus providing a tool for studying cosmology. We find good agreement between the mass map and 
the distribution of massive galaxy clusters identified using a red-sequence cluster finder. Potential candidates 
for super-clusters and voids are identified using these maps. We measure the cross-correlation between the 
mass map and a magnitude-limited foreground galaxy sample and find a detection at the 6.8 <t level with 20 
arcminute smoothing. These measurements are consistent with simulated galaxy catalogs based on ACDM N- 
body simulations, suggesting low systematics uncertainties in the map. We summarize our key findings in this 
letter; the detailed methodology and tests for systematics are presented in a companion paper. 
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PACS numbers: 


INTRODUCTION 

Gravitational lensing refers to the bending of light due to 
the curvature of space-time induced by massive bodies [ 1 ]. 
This effect allows one to probe the total matter distribution in 
the Universe, including both luminous and dark matter. Weak 
lensing is the technique of using the subtle gravitational lens¬ 
ing effect of a large number of galaxies to statistically infer 
the large-scale matter distribution in the Universe [see 2, 3, 
for detailed reviews]. The measurement is based on small, 
percent-level “shears”, or distortions of galaxy shapes due to 
lensing. With several ongoing large optical surveys collecting 
data [4-7], this technique is one of the most powerful probes 
for constraining the nature of dark energy [ 8 ]. 

Conventional weak lensing analyses involve calculating the 
N-point statistics of the shear field. In particular, the cosmic 
shear measurement, which refers to the 2 -point correlation 
function of the shear field in configuration space, has been 
measured in several earlier datasets [9-14]. Shear 7 is de¬ 
fined to be a combination of second derivatives of the lensing 
potential \j/. 


7 = 7i +'72 = ^ iV.ii - ¥,22)+i¥,i2, (1) 

“¥,ij = ^^¥ is the second partial derivative with re¬ 

spect to the angular sky coordinates 0 , of \j/ (assuming a spa¬ 
tially flat Universe in the Newtonian limit of GR), which is 
defined as [15] 

Xlf{e,r)=2 (2) 

Jo rr ' ' 

In the above equation, r is the comoving distance and <I> is the 
3D gravitational potential, whose spatial structure and time 
evolution contains cosmological information. 

Instead of measuring statistics based on shear, here we fo¬ 
cus on an alternative approach by converting shear into the 
projected density field, the convergence K, also a combination 
of second derivatives of \j/, 

K=^¥^¥=\{¥,n + ¥,22)- (3) 

The convergence directly represents the integrated mass distri¬ 
bution, which can be seen by using the cosmological Poisson 
equation and the Limber approximation to re-write Eq. (3) as 
[ 2 ] 


K{e,r) 


3H^Q.m r ,r'{r-r') d {6, r') 

2 Jo r fl(r') 


(4) 


where Hq is the Hubble constant today, is the total matter 
density today, a is the cosmological scale factor, and 5 = (A — 
A)/A is the mass overdensity (A and A are the 3D density and 
mean density respectively). In practice we integrate over the 


redshift distribution of source galaxies as shown in Eq. 15 of 
the accompanying paper [16]. 

Note that the same weak lensing effect also introduces 
distortions in the observed cosmic microwave background 
(CMB) maps. Reconstructing the convergence map from the 
CMB gives the integrated mass up to the surface of last scat¬ 
tering (z ~ 1100). Compared to the weak lensing convergence 
map constructed from galaxies, the CMB convergence map 
typically covers a larger area with lower spatial resolution, 
and the sources of the lensing effect (the CMB photons) come 
from a single redshift plane [17-19]. In this letter, we use 
“weak lensing mass maps” to refer to convergence maps gen¬ 
erated from source galaxies. 

Weak lensing mass maps supplement measurements based 
on shear in many ways. Mass maps can be easily cross- 
correlated with other data since they represent a scalar, the 
local (projected) mass density, while the shear is a complex 
variable and is sensitive to the global mass distribution. Cross 
correlating with X-ray and Sunyaev-Zel’dovich observations 
helps us understand the relation of hot gas and dark matter 
in galaxy clusters. Cross correlating with the CMB conver¬ 
gence map provides an important cross check of lensing mea¬ 
surements using different tracers. Other applications of mass 
maps include peak statistics [20-24], higher-order moments 
of K [25], and the identification of superclusters and cosmic 
voids [26]. 

The methodology of generating weak lensing mass maps 
has been demonstrated in earlier work. Massey et al. [27] 
generated a 3D mass map using COSMOS data in a 1.64 deg^ 
area. The high-quality shear measurements and redshift in¬ 
formation allow for good mass reconstruction on small scales 
and in the radial direction. Van Waerbeke et al. [28], on the 
other hand, focused on larger-scale information and gener¬ 
ated 2D wide-held mass maps from four helds of size 25-72 
deg^ in the Canada-Erance-Hawaii Telescope Lensing Survey 
(CEHTLenS). Our work is similar to Van Waerbeke et al. [28], 
but uses one contiguous region of 139 deg^ from the Dark En¬ 
ergy Survey [DES, 5, 29] data. This is the hrst step towards 
building mass maps from the full DES data set. 

The data used in this work is part of the Science Verihca- 
tion (SV) dataset from DES, an ongoing ground-based galaxy 
survey that is scheduled to operate from September 2013 to 
Eebruary 2018. The SV data were collected between Novem¬ 
ber 2012 and Eebruary 2013 shortly after the commissioning 
of the new wide-held mosaic camera, the Dark Energy Cam¬ 
era [DECam, 30-32] on the 4m Blanco telescope at the Cerro 
Tololo Inter-American Observatory (CTIO) in Chile. This 
data was used to test survey operations and assess data qual¬ 
ity. The images are taken in 5 optical hlter bands (grizY) on 
a total area of ~ 250 deg^ and reach close to the expected full 
depth of DES at 23.9. 

The main goal of this work is to reconstruct the weak 
lensing mass map from shear measurements of the DES SV 
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TABLE I: Catalogs and selection criteria used to construct the fore¬ 
ground and background sample for this work, and the number of 
galaxies in each sample after all the cuts are applied. The redshift 
cut is based on the mean redshift output from the BPZ photo-z code 
and the magnitude cut is based on the MAG_AUTQ parameter in the 
S Extractor output. 



Background (source) 

Foreground (lens) 

Input catalog 

ngmixOll imSshape 

SVAl Gold 

Photo-z 

0.6<z<L2 

0.1<z<0.5 

Selection 

“conservative additive” 

i<22 

Number of galaxies 

1,111,487 1,013,317 

1,106,189 

Number density 
(arcmin^^) 

2.22 

2.03 

2.21 

Mean redshift 

0.826 

0.825 

0.367 


data in a 139 deg^ contiguous region overlapping with the 
South Pole Telescope Survey (the SPT-E field). We present 
the methodology used for the map construction, followed by 
cross-correlation results and conclusions. Throughout the pa¬ 
per, we adopt the following cosmological parameters: = 

0.3, Oa = 0.7, = 0.0, h — 0.72. A detailed account of this 

work can be found in a companion paper in PRD [16]. 


METHODOLOGY 
Data and simulations 

Our galaxy samples are based on the DES SV Gold cata¬ 
log (Rykoff et ak, in preparation) and several extensions to 
it. The Gold catalog is a product of the DES Data Man¬ 
agement [DESDM, 33-36] pipeline version “SVAl” (Yanny 
et ak, in preparation), which includes calibrated photometry 
and astrometry, object morphology, object classification and 
masking of the co-add SV images. DESDM utilizes the soft¬ 
ware packages SCAMP [37], SWARP [38], PSEEx [39] and 
SExtractor [40] in the pipeline. 

Several additional catalogs are used in this work. We use a 
photometric redshift (photo-z) catalog from the photo-z code 
Bayesian Photometric Redshifts [BPZ, 41, 42]. We use two 
shear catalogs from the ngmix code [43] and the imSshape 
code [44]. The two independent shear catalogs allows us to 
assess the robustness of the measurement. The shear mea¬ 
surement algorithms operate on single-exposure images and 
measure the galaxy shapes, or “ellipticities”, by jointly fitting 
the images of the same galaxies obtained in different expo¬ 
sures with one galaxy model and the different point-spread- 
function (PSP) model in each image. The resulting ellipticity 
is a noisy estimator for shear [2]. The shear estimates used in 
this work have been tested rigorously as described in Jarvis et 
al. (in preparation). 

We extract from these catalogs background (“source”) and 
foreground (“lens”) galaxy samples. The objective is to con¬ 
struct the convergence, or mass map, from the background 
sample and cross-correlate it with the weighted galaxy map 
built from the foreground sample. Table I lists the final se¬ 


lection criteria for the samples. The foreground sample is 
magnitude-limited at i = 22 , while the background sample 
is selected through a series of lensing tests (Jarvis et al. in 
preparation) and is not complete. The incompleteness of the 
background sample affects only the spatial distribution of the 
noise on these maps but does not bias the signal. In the com¬ 
panion paper we describe in detail the construction of these 
samples and also discuss a second foreground sample com¬ 
posed of luminous red galaxies (LRGs). Note that the plots 
in this letter rely on the ngmix shear catalog. However we 
analyzed both shear catalogs to assess their statistical consis¬ 
tency. The “conservative additive” selection criteria on the 
background sample involves a combination of signal-to-noise 
(S/N) cuts, size cuts and other quality cuts. 

To facilitate our understanding of possible systematics in 
the procedure of constructing the mass map, we use a set 
of simulated galaxy catalogs that we match closely to the 
characteristics of the data (including intrinsic galaxy prop¬ 
erties, galaxy number counts, noise, photo-z errors, survey 
mask). We use the simulated galaxy catalogs developed for 
the DES collaboration [45]. The catalog is based on three flat 
ACDM dark matter-only N-body simulations with different 
resolutions. Galaxies are populated using the prescriptions de¬ 
rived from a high-resolution simulation using SubHalo Abun¬ 
dance Matching techniques [45^7]. Photometric properties 
for each galaxy are then assigned so that the magnitude-color- 
redshift distribution reproduces that observed in the SDSS 
DR 8 [48] and DEEP2 [49] data. Weak lensing parameters 
(shear and convergence) are assigned to each galaxy based on 
the high-resolution ray-tracing algorithm Curved-sky grAvi- 
tational Lensing for Cosmological Light conE simulatioNS 
[50]. Details of the data and simulation catalogs are presented 
in the companion paper. 


Mass and weighted galaxy maps 


Eq. (1) and Eq. (3) can be Eourier transformed to get a sim¬ 
ple relationship between the Eourier transforms of the shear 
and convergence, denoted 7 and k [51 ]: 


f j — ^2 +2!Tif2 


(5) 

( 6 ) 


where are the components of the angular wavenumber. The 
above equations hold for 1>Q. 

In practice, we pixelate the shear measurements into a map 
of 5x5 arcmin^ pixels and Eourier transform the map. We 
then use Eq. (5) to obtain k and inverse Eourier transform to 
yield our final real-space convergence map. In an ideal sce¬ 
nario, this reconstructed convergence map does not contain 
an imaginary component. However, due to noise, the finite 
area of the map, and masking, a non-zero imaginary compo¬ 
nent is recovered. We separate the real and imaginary parts 
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FIG. 1: The DES SV weighted foreground galaxy maps Kg ,„ain (left). E-mode convergence map Ke (middle) and B-mode convergence map 
Kb (right) are shown in these panels. All maps are generated with 5x5 arcmin^ pixels and 20 arcmin RMS Gaussian smoothing. In the Kg and 
Ke maps, red areas corresponds to overdensities and blue areas to underdensities. White regions correspond to the survey mask. The scale of 
the Gaussian smoothing kernel is indicated by the Gaussian profile on the upper right comer of the right panel. The Ke map is overlaid by 
Redmapper galaxy clusters with optical richness A > 20. The radius of the circles scale with A. The black and white squares show the super 
cluster and super void candidate we investigate in Fig. (2). 


of the measured convergence map into E- and B-modes, or 
K — Ke + iKb- The B-mode convergence is a useful diagnostic 
tool for testing systematics, as it should vanish for real lensing 
signals on a sufficiently large area. Finally, as the uncertainty 
in this reconstruction is formally infinite for a discrete set of 
noisy shear estimates, it is important to apply a filter to re¬ 
move the high-frequency noise [52]. In this work we apply a 
Gaussian filter of different sizes. In the companion paper [16] 
we use simulations to quantify the degradation in Ke and the 
level of Kb expected from the noise and masking in the data. 
We find that our results are consistent with that expected from 
simulations. 


One of the main goals of this work is to cross-corTelate the 
mass map with the foreground galaxy distribution. For this 
purpose, we construct a second mass map assuming that the 
foreground galaxy sample traces the mass distribution - we re¬ 
fer to this map as Kg. It is constructed using equation Eq. (4) 
with 5 replaced by 5g, the fractional overdensity of galaxy 
counts. Under the assumption of linear bias (i.e. galaxy over¬ 
densities are linearly proportional to the total mass overdensi¬ 
ties, which is expected to be valid on sufficiently large scales), 
the smoothed Kg is simply a product of the mass map K with 
a constant bias factor. For our foreground galaxy sample, the 
linear bias is valid above 5-10 arcmin scales, which is the fo¬ 
cus of our study [53]. In practice, the limited redshift range of 
our foreground galaxy sample means that we cannot expect a 
perfect estimate of the mass map even if the bias factor were 
unity. 


RESULTS 

Fig. (1) shows the resulting weighted galaxy map and the E 
and B-mode convergence maps generated from the procedure 
described above. The maps shown are for a Gaussian smooth¬ 
ing of 20 arcmin RMS. We expect Ke to correlate with Kg, 
while Kb should not correlate with either of the other maps. 

Correlation with clusters 

The Ke map shown in the middle panel of Fig. (1) is over¬ 
laid with galaxy clusters detected in the same data using the al¬ 
gorithm Redmapper [54]. Each cluster is represented by a cir¬ 
cle with radius proportional to the optical richness A, which is 
related to mass via a roughly linear relation (see Rykoff et al. 
[54] for details of the mass calibration of A). We select only 
clusters with A > 20, which corresponds to mass > 1.7 x 10*^ 
M© (A = 80 corresponds to mass ~ 7.6 x lO*"'' M©). Visu¬ 
ally, one can see that the spatial distribution of the clusters 
traces the mass map very well, with most clusters detected in 
or around the high Ke regions. 

We analyze the redshift distributions of the clusters in the 
high and low mass density regions. Two examples are shown 
in Fig. (2), where we plot (in blue) the lensing efficiency 
and A-weighted redshift distribution of the clusters within a 
1 degree radius of the identified high and low-mass positions 
marked in Fig. (1). Compared to the average redshift distribu¬ 
tion of clusters (overlaid in grey), we find that the high-mass 
(low-mass) regions indeed contain many more (fewer) clusters 
than average. The redshift binning is Az=0.03, corresponding 
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FIG. 2: Blue lines show the richness-weighted redshift distribution 
of redmapper galaxy clusters along along overdense (left) and under- 
dense (right) regions in the convergence map marked by the black 
and white squares in Fig. (1). The (RA, Dec) positions of the each 
region is shown in the upper right comer of each panel. The thick 
grey line shows the average redshift distribution over the full map. 
Both lines are weighted by the lensing efficiency. The redshift range 
above z = 0.6 (marked with the shaded grey area) overlap with the 
background sample, hence the interpretation of the structures there is 
more complicated. 


to between 1.5-3 CJ^ in this redshift range, where O). is the clus¬ 
ter photo-z error uncertainty. The photo-z’s for Redmapper 
clusters are very well determined (ff^ « 0.01(1 +z)), which 
is important for the identihcations of the 3D structures. Us¬ 
ing these histograms we can identify potential candidates for 
super-clusters. For example, the peak at z ~ 0.14 in the left 
panel parked in red indicates that this spatial structure is con¬ 
tained in a redshift range localized to within about 100 Mpc 
along the line of sight. This line of sight has multiple struc¬ 
tures at different redshifts, others have just one or two. The 
redshift range above z = 0.6 (marked with the shaded grey 
area) overlaps with the background sample, hence the inter¬ 
pretation of their relation with the mass map is more com¬ 
plicated. The largest mass concentrations are investigated in 
more detail in the companion paper and in follow-up studies. 


Mass-galaxy correlation 

Next, we investigate quantitatively the correlation between 
the foreground galaxies and the mass map by calculating the 
Pearson correlation coefficient between the two maps over a 
range of smoothing scales that span 5 to 40 arcmin. That is, 
we calculate 




<yKE<^Kg 


(7) 


where {KEKg) is the covariance between Ke and Kg, and 0 ,^^ 
and are the standard deviations of the two maps. In this 
calculation, pixels in the masked region are not used. We also 
remove pixels within 10 arcmin of the boundaries to avoid 
signihcant artefacts from the smoothing. Similarly we cal¬ 
culate the Pearson correlation coefficient between Kg and the 
other maps to check for any signihcant systematic effects. The 
errors on the correlation coefficients are estimated by a jack¬ 



FIG. 3: The Pearson correlation coefficient between the foreground 
galaxy and convergence maps is shown as a function of smoothing 
scale. The solid and open symbols show correlation coefficients from 
the E and B-modes of the convergence respectively. The grey shaded 
regions show the IcT bounds from simulations for the correlation 
between the E and B-mode convergence and the foreground galax¬ 
ies, with the same pixelization and smoothing as the data as well as 
sources of statistical uncertainty. The green points show the correla¬ 
tion between E and B-modes of the convergence map. The various 
correlation coefficients with the B-mode convergence are consistent 
with zero. Uncertainties on all measurements are estimated using 
jackknife resampling. 


knife resampling of 10 deg^ sub-regions of the maps (each 
jackknife subsample is ^ 93% of the total area). 

The results are shown in Fig. (3). We find that the Pear¬ 
son correlation coefficient between Kg and the E-mode conver¬ 
gence is 0.39 ±0.06 at 10 arcmin smoothing and 0.52 ±0.08 
at 20 arcmin smoothing. This corresponds to a ^ 6.8(7 sig¬ 
nificance at these scales. The correlation between the B-mode 
convergence and the Kg maps is consistent with zero at all 
smoothing scales. The correlation between the E and B-modes 
convergence is also consistent with zero. The grey shaded re¬ 
gions show the 1 a range of results from the simulated galaxy 
catalogs modelled to match the main characteristics of the data 
samples. The black data points agree well with the simula¬ 
tions, suggesting there are no significant contributions to our 
signal from systematic errors. 

To further examine the potential contamination by system- 
atics in the maps, we construct maps of 20 quantities associ¬ 
ated with the observing conditions (e.g. airmass, extinction, 
seeing, PSE ellipticity etc.) and cross correlate with our Ke 
and Kg maps. We find that none of these quantities contribute 
significantly to the cross correlation signal we have measured, 
with most of them consistent with zero. Details are presented 
in the companion paper. 


SUMMARY 


We present in this letter a weak lensing convergence map 
generated from shear measurements in the 139 deg^ SPT-E 
field in the Dark Energy Survey Science Verification data. The 
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mean redshift of the source galaxies is 0.82 and corresponds 
to a comoving distance of 2.9 Gpc. This map probes the pro¬ 
jected total mass (luminous and dark), with matter approxi¬ 
mately half-way between us and the source galaxies making 
the most contribution to the lensing. We study the correla¬ 
tion of the mass map with galaxies and clusters that trace the 
foreground mass distribution. 

The spatial distribution of galaxy clusters identified in the 
same data using an independent technique is highly correlated 
with the mass map. The combination of the mass map and 
the cluster catalog provide a powerful tool for exploring po¬ 
tential super-clusters and super-voids in the Universe. Cross- 
correlating the E-mode mass map with a magnitude-limited 
foreground galaxy sample gives a 6.8(7 detection at 20 ar- 
cminute smoothing, while the cross correlation between B- 
mode mass map and the galaxies is consistent with zero on 
all scales. The cross-correlation between E and B-mode mass 
map are also consistent with zero. These results are consistent 
with simulations of the ACDM model in which we have mod¬ 
eled several sources of statistical uncertainties in the lensing 
and weighted galaxy maps. More detailed analysis, simula¬ 
tion and systematics tests are described in a companion PRD 
paper [16]. 

Topics for follow-up studies include the study of galaxy 
bias, identification of super-clusters and super-voids, higher 
order moments of the mass map, and cross-correlation with 
the CMB and other observations. With the full set of data 
from DES in a few years (~ 35 times the size of the SV data 
used in this work), we expect the mass maps to be a powerful 
tool for cosmology. 
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